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Abstract 

We present distribution independent bounds on the generalization misclassification perfor- 
mance of a family of kernel classifiers with margin. Support Vector Machine classifiers 
(SVM) stem out of this class of machines. The bounds are derived through computations 
of the Ky dimension of a family of loss functions where the SVM one belongs to. Bounds 
that use functions of margin distributions (i.e. functions of the slack variables of SVM) are 
derived. 
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1 Introduction 

Deriving bounds on the generalization performance of kernel classifiers has been an important 
theoretical topic of research in recent years [4, 8, 9, 10, 12]. We present new bounds on the gen- 
eralization performance of a family of kernel classifiers with margin, from which Support Vector 
Machines (SVM) can be derived. The bounds use the V^ dimension of a class of loss functions, 
where the SVM one belongs to, and functions of the margin distribution of the machines (i.e. 
functions of the slack variables of SVM - see below). 
We consider classification machines of the form: 

min EZiV(yiJ(xi)) 

subject to \\f\\ 2 K <A 2 (1) 

where we use the following notation: 

• D m = {(xi,?/i), . . . , (x m ,y m )}, with (x^y*) e R n x {-1, 1} sampled according to an un- 
known probability distribution P(x,y), is the training set. 

• V(y: /( x )) is the loss function measuring the distance (error) between /(x) and y. 

• f is a function in a Reproducing Kernel Hilbert Space (RKHS) 7i defined by kernel K, 
with ||/ 1| |f being the norm of / in Ti [11, 2]. We also call / a hyperplane, since it is such 
in the feature space induced by the kernel K [11, 10]. 

• A is a constant. 

Classification of a new test point x is always done by simply considering the sign of /(x). 
Machines of this form have been motivated in the framework of statistical learning theory. We 
refer the reader to [10, 6, 3] for more details. In this paper we study the generalization perfor- 
mance of these machines for choices of the loss function V that are relevant for classification. In 
particular we consider the following loss functions: 

• Misclassification loss function: 

V(y, /(x)) = V msc (yf(x)) = 0(-y/(x)) (2) 

• Hard margin loss function: 

V(y, /(x)) = V hm (yf(x)) = 6(1 - y/(x)) (3) 

• Soft margin loss function: 

V(y, /(x)) = V sm (yf(s)) = 6(1 - y/(x))(l - y/(x)), (4) 

where 6 is the Heavyside function. Loss functions (3) and (4) are "margin" ones because the only 
case they do not penalize a point (x, y) is if yf(x) > 1. For a given /, these are the points that 
are correctly classified and have distance 4^r > ^fht from the surface /(x) = (hyperplane in 
the feature space induced by the kernel K [10]). For a point (x, y), quantity Hffl- is its margin, 




Figure 1: Hard margin loss (line with diamond-shaped points), soft margin loss (solid line), 
nonlinear soft margin with a = 2 (line with crosses), and a — | (dotted line) 

and the probability of having Hjrr > S is called the margin distribution of hypothesis /. For 

SVM, quantity 0(1 — yj/(xj))(l — t/j/(xj)) is known as the s/ac/c variable corresponding to training 

point (x^) [10]. 

We will also consider the following family of margin loss functions (nonlinear soft margin loss 

functions): 

V(y, /(x)) = V*(y/(x)) = 0(1 - y/(x))(l - y/(x))*. (5) 

Loss functions (3) and (4) correspond to the choice of a — 0, 1 respectively. In figure 1 we plot 
some of the possible loss functions for different choices of the parameter a. 



To study the statistical properties of machines (1) we use some well known results that we now 
briefly present. First we define some more notation, and then state the results from the literature 
that we will use in the next section. 
We use the following notation: 

• R% mp (f) = YhLi V{yi-> /( x i)) i s the empirical error made by / on the training set D m , using 
V as the loss function. 

• R v (f) = Jjf. x ;_in V"(y, /(x)) P(x, y) dx. dy is the expected error of / using V as the loss 
function. 



• Given a hypothesis space of functions T (i.e. T = {/ G 7i : ||/|| 2 < A 2 }), we note by hY T 
the V^ dimension of the loss function V(y, /(x)) in J 7 , which is defined as follows [1]: 

Definition 1.1 Let A < V(y, /(x)) < B, f e T , with A and B < oo. The V 1 - dimension of V 
in JF (of the set of functions {V(y, /(x)) | / G -F}j zs defined as the the maximum number h 
of vectors (xi, yi) . . . , (x^, y/j) i/iat can be separated into two classes in all 2 h possible ways using 
rules: 



class 1 if: V(y t , /(x^)) > s + 7 
c/ass -1 */• V(j/i, /(xj)) < s - 7 

/or / G JF and some s > 0. 7/j /or any number m, it is possible to find m points (x 1; yi) . . . , (x m , y m ] 
that can be separated in all the 2 m possible ways, we will say that the V^-dimension ofV in T is 
infinite. 

If instead of a fixed s for all points we use a different Si for each (x,, y^), we get what is called the 
fat-shattering dimension fat 7 [1]. Notice that definition (1.1) includes the special case in which 
we directly measure the V^ dimension of the space of functions F, i.e. V(y, /(x)) = /(x). We 
will need such a quantity in theorem 2.2 below. 

Using the V 1 dimension we can study the statistical properties of machines of the form (1) based 
on a standard theorem that characterizes the generalization performance of these machines. 

Theorem 1.1 (Alon et al., 1993) Let A < V(y,f(x)) < B , f E T ', T be a set of bounded 
functions. For any e > 0, for all m > -^ we have that if hX T is the V^ dimension of V in T for 
7 = ae (a > -^), h^ r finite, then: 



Pr jsup \R^ mp (f) - R v (f)\ >e\< Q{e, m, h?) 



(6) 



where Q is an increasing function of hY T and a decreasing function of e and m, with Q — > as 



7 
m — > 00. 



In [1] the fat-shattering dimension was used, but a close relation between that and the V^ 
dimension [1] make the two equivalent for our purpose 1 . Closed forms of Q can be derived (see 
for example [1]) but we do not present them here for simplicity of notation. Notice that since 
we are interested in classification, we only consider e < 1, so we will only discuss the case 7 < 1 
(since 7 is about ^e). 

In "standard" statistical learning theory the VC dimension is used instead of the V 1 one [10]. 
However, for the type of machines we are interested in the VC dimension turns out not to be 
appropriate: it is not influenced by the choice of the hypothesis space T through the choice of A, 
and in the case that T is an infinite dimensional RKHS, the VC-dimension of the loss functions 
we consider turns out to be infinite (see for example [5]). Instead, scale-sensitive dimensions 
(such as the V^ or fat-shattering one [1]) have been used in the literature, as we will discuss in 
the last section. 

2 Main results 

We study the loss functions (2-5). For classification machines the quantity we are interested 
in is the expected misclassification error of the solution / of problem 1. With some abuse of 
notation we note this with R msc . Similarly we will note with R hm , R sm , and R a the expected 
risks using loss functions (3), (4) and (5), respectively, and with R^ p , R 8 ^, and R° mp , the 



^n [1] it is shown that V 1 < fat 7 < ±Vi 



corresponding empirical errors. We will not consider machines of type (1) with V msc as the loss 
function, for a clear reason: the solution of the optimization problem: 

min TT=iK-Vif^i)) 
subject to \\f\\ 2 K <A 2 

is independent of A, since for any solution / we can always rescale / and have the same cost 

E^-i/i/CxO). 

For machines of type (1) that use V sm or V a as the loss function, we prove the following: 

Theorem 2.1 The V^ dimension h for 6(1 — yf(x))(l — yf (x)) CT in hypothesis spaces J- a = {/ G 
H\\\f\\ 2 K < A 2 } (of the set of function {9(l-yf(x))(l-yf(x))« | / e T A }) and y e {-1, 1}, 

is finite for V < 7. If D is the dimensionality of the RKHS 7i, R 2 is the radius of the smallest 
sphere centered at the origin containing the data x in the RKHS, and B > 1 is an upper bound 
on the values of the loss function, then h is upper bounded by: 



0(min(D, !££-)) for a < 1 
U(mm(V, i )) for a > 1 



Proof 

The proof is based on the following theorem [7] (proved for the fat-shattering dimension, but as 
mentioned above, we use it for the "equivalent" V^ one). 

Theorem 2.2 [Gurvits, 1997] The V^ dimension h of the set of functions 2 Ta = {/ € T^IH/Hl; < 
A 2 } is finite for V 7 > 0. If D is the dimensionality of the RKHS, then h < 0(min(D, ^-§-)), 
where R 2 is the radius of the smallest sphere in the RKHS centered at the origin here the data 
belong to. 

Let 2N be the largest number of points {(xi, yi), . . . , (x2jv, V2n)} that can be shattered using the 

rules: 

class 1 if $(l-y i f(x i ))(l-y i f(x i )y>s + 1 

class -lif e(l-y t f(^))(l-y t f(^)r<s--f { ' > 

for some s with < 7 < s. After some simple algebra these rules can be decomposed as: 

class 1 if /(xj) - 1 < -(s + 7)- (for y { = 1 ) 

or /(xi) + 1 > (s + 7) a (for Vi = -1 ) 

class - 1 if /(xj) - 1 > -(s - 7)- (for ^ = 1 ) 

or /(xi) + 1 < (s - 7)- (for j/ 4 = -1 ) 

From the 2N points at least N are either all class -1, or all class 1. Consider the first case (the 
other case is exactly the same), and for simplicity of notation let's assume the first N points are 



2 As mentioned above, in this case we can consider V(y, /(x)) = /(x). 



class -1. Since we can shatter the 2N points, we can also shatter the first N points. Substituting 
Hi with 1, we get that we can shatter the N points {x 1; . . . , x^} using rules: 

class 1 if /(xi) + l > (s + 7)- , g x 

class — 1 if /(xj) + 1 < (s — 7)" 

Notice that the function /(xj) + 1 has RKHS norm bounded by A 2 plus a constant C (equal to 
the inverse of the eigenvalue corresponding to the constant basis function in the RKHS - if the 
RKHS does not include the constant functions, we can define a new RKHS with the constant 
and use the new RKHS norm). Furthermore there is a "margin" between (s + 7) ^ and (s — 7)^ 
which we can lower bound as follows. 
For a < 1, assuming - is an integer (if not, we can take the closest lower integer), 

I(( s + 7 )i- ( a _ 7 )*) = I(( s + 7 )_( a _ 7 )) j£( s + 7 )*-i-*( s _ 7 )*] > 77 -i = 7 -. ( io) 

For a > 1, a integer (if not, we can take the closest upper integer) we have that: 

2 7 = ((* + #) CT - - #) CT = ((* + # - (a - #) (l> + T)-)"- 1 -*^ - #) fc ) < 

\fc=0 / 

< {{s + 7) ct — (s — r ))^)aB~^~ 
from which we obtain: 

^ (> + #-(*-#)>— kr (11) 

Therefore N cannot be larger than the V^ dimension of the set of functions with RKHS norm 
< A 2 + C and margin at least 7" for a < 1 (from eq. (10)) and — ^— T for <r > 1 (from eq. (11)). 

Using theorem 2.2, and ignoring constant factors (also ones because of C), the theorem is proved. 

□ 



In figure 2 we plot the V^ dimension for R 2 A 2 = 1, B = 1, 7 = 0.9, and D infinite. Notice 
that as a — > 0, the dimension goes to infinity. For a = the V^ dimension becomes the same 
as the VC dimension of hyperplanes, which is infinite in this case. For a increasing above 1, the 
dimension also increases: intuitively the margin 7 becomes smaller relatively to the values of the 
loss function. 
Using theorems 2.1 and 1.1 we can bound the expected error of the solution / of machines (1): 

Pr{\RL P (f) ~ R V (f)\ >*}< Sbm, hj, (12) 

where V is V sm or V a . To get a bound on the expected misclassification error R msc (f) we use 
the following simple observation: 

U ms %,/(x)) < U%,/(x)) for Va, (13) 




Sigma 



Figure 2: Plot of the V^ dimension as a function of a for 7 = .9 

So we can bound the expected misclassification error of the solution of machine (1) under V sm 
and V 7 using the V n dimension of these loss functions and the empirical error of / measured using 
again these loss functions. In particular we get that for V<7, with probability 1 — Q(e, m, hj r ): 

R mS V)<R a em P (f)+e (14) 

where e and 7 are related as stated in theorem 1.1. 



Unfortunately we cannot use theorems 2.1 and 1.1 for the V hm loss function. For this loss 
function, since it is a binary-valued function, the V^ dimension is the same as the VC-dimension, 
which, as mentioned above, is not appropriate to use in our case. Notice, however, that for 
a — > 0, V a approaches V pointwise (from theorem 2.1 the V^ dimension also increases towards 
infinity). Regarding the empirical error, this implies that R a — > R hm , so, theoretically, we can 

using: 



still bound the misclassification error of the solution of machines with V 



R" 



V) < ReZ(f) 



e + max(i^ (f) 



RJ2 v (f),0), 



emp\ 



(15) 



where R° mp (f) is measured using V a for some a. Notice that changing a we get a family of 
bounds on the expected misclassification error. Finally, we remark that it could be interesting to 
extend theorem 2.1 to loss functions of the form 0(1 —yf(x))h(l —yf(x)), with h any continuous 
monotone function. 



3 Discussion 

In recent years there has been significant work on bounding the generalization performance of 
classifiers using scale-sensitive dimensions of real-valued functions out of which indicator func- 
tions can be generated through thresholding (see [4, 9, 8], [3] and references therein). This is 
unlike the "standard" statistical learning theory approach where classification is typically stud- 
ied using the theory of indicator functions (binary valued functions) and their VC-dimension [10]. 



The work presented in this paper is similar in spirit with that of [3], but significantly different 
as we now briefly discuss. 

In [3] a theory was developed to justify machines with "margin". The idea was that a "better" 
bound on the generalization error of a classifier can be derived by excluding training examples 
on which the hypothesis found takes a value close to zero (as mentioned above, classification 
is performed after thresholding a real valued function). Instead of measuring the empirical 
misclassification error, as suggested by the standard statistical learning theory, what was used 
was the number of misclassified training points plus the number of training points on which the 
hypothesis takes a value close to zero. Only points classified correctly with some "margin" are 
considered correct. In [3] a different notation was used: the parameter A in equation (1) was 
fixed to 1, while a margin tp was introduced inside the hard margin loss, i.e 6(ip — yf(x)). Notice 
that the two notations are equivalent: given a value A in our notation we have ip = A~ l in the 
notation of [3]. Below we adapt the results in [3] to the setup of this paper, that is, we set ip — 1 
and let A vary. Two main theorems were proven in [3]. 

Theorem 3.1 (Bartlett, 1998) For a given A, with probability 1 — S, every function f with 
\ 2 K < A 2 has expected misclassification error R msc (f) bounded as: 



Ti ■(/) < Ti'ZW - \l^(dln(Mem/d) log 2 (578m) + ln(4/5), (16) 



where d is the fat-shattering dimension f at^ of the hypothesis space {/ : ||/||#- < A } for 7 = -^. 

Unlike in this paper, in [3] this theorem was proved without using theorem 1.1. Although prac- 
tically both bound (16) and the bounds derived above are not tight and therefore not practical, 
bound (16) seems easier to use than the ones presented in this paper. 

It is important to notice that, like bounds (12), (14), and (15), theorem 3.1 holds for a fixed A 
[3]. In [3] theorem 3.1 was extended to the case where the parameter A (or ip in the notations of 
[3]) is not fixed, which means that the bound holds for all functions in the RKHS. In particular 
the following theorem gives a bound on the expected misclassification error of a machine that 
holds uniformly over all functions: 

Theorem 3.2 (Bartlett, 1998) For any f with WfWx < oo 7 with probability 1 — S, the mis- 
classification error R mcs (f) of f is bounded as: 



R msc (f) < R^ p (f) + J-(dln(34em/d) log 2 (578m) + ln(8\\f\\/5), (17) 

where d is the fat-shattering dimension fat^ of the hypothesis space consisting of all functions in 
the RKHS with norm < ||/|||-, and with 7 = 32 n ,,, . 

Notice that the only differences between (16) and (17) are the Zn(8||/||/5) instead of ln(4/S), and 
that 7 = 32 n ,n instead of 7 = -^j. 

So far we studied machines of the form (1), where A is fixed a priori. In practice learning 
machines used, like SVM, do not have A fixed a priori. For example in the case of SVM the 
problem is formulated [10] as minimizing: 

min T£ieQ--Vif(xi))0--Vif(*i)) + M\f\\ 2 K (18) 



where A is known as the regularization parameter. In the case of machines (18) we do not know 
the norm of the solution \\f\\ 2 K before actually solving the optimization problem, so it is not clear 
what the "effective" A is. Since we do not have a fixed upper bound on the norm \\f\\\ a priori, 
we cannot use the bounds of section 2 or theorem 3.1 for machines of the form (18). Instead, 
we need to use bounds that hold uniformly for all A (or ip if we follow the setup of [3]), for 
example the bound of theorem 3.2, so that the bound also holds for the solution of (18) we find. 
In fact theorem 3.2 has been used directly to get bounds on the performance of SVM [4]. A 
straightforward applications of the methods used to extend theorem 3.1 to 3.2 can also be used 
to extend the bounds of section 2 to the case where A is not fixed (and therefore hold for all / 
with ll/H < oo), and we leave this as an exercise. 

There is another way to see the similarity between machines (1) and (18). Notice that the 
formulation (1) the regularization parameter A of (18) can be seen as the Lagrange multiplier 
used to solve the constrained optimization problem (1). That is, problem (1) is equivalent to: 

m 

max A mi n/ JTVfo, /(*)) + X(\\f\\ 2 K - A 2 ) (19) 

8=1 

for A > 0, which is similar to problem (18) that is solved in practice. However in the case of 
(19) the Lagrange multiplier A is not known before having the training data, unlike in the case 
of (18). 

So, to summarize, for the machines (1) studied in this paper, A is fixed a priori and the "regular- 
ization parameter" A is not known a priori, while for machines (18) the parameter A is known a 
priori, but the norm of the solution (or the effective A) is not known a priori. As a consequence 
we can use the theorems of this paper for machines (1) but not for (18). To do the second we 
need a technical extension of the results of section 2 similar to the extension of theorem 3.1 to 
3.2 done in [3]. On the practical side, the important issue for both machines (1) and (18) is 
how to choose A or A. We believe that the theorems and bounds discussed in sections 2 and 
3 cannot be practically used for this purpose. Criteria for the choice of the regularization pa- 
rameter exist in the literature - such as cross validation and generalized cross validation - (for 
example see [10, 11], [6] and references therein), and is the topic of ongoing research. Finally, 
as our results indicate, the generalization performance of the learning machines can be bounded 
using any function of the slack variables and therefore of the margin distribution. Is it, however, 
the case that the slack variables (margin distributions or any functions of these) are the quan- 
tities that control the generalization performance of the machines, or there are other important 
geometric quantities involved? Our results suggest that there are many quantities related to the 
generalization performance of the machines, but it is not clear that these are the most important 
ones. 
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